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' We show that the Hidden Local Symmetry Model, supplemented with 

well-known procedures for breaking flavor SU(3) and nonet symmetry, pro- 
vides all the information contained in the standard Chiral Perturbation The- 
, ory (ChPT) Lagrangian C^^^ + C^^^ . This allows to rely on radiative decays 

Qj [ of light mesons {VP^ and -P77) in order to extract some numerical informa- 

' tion of relevance to ChPT: a value for Ai = 0.20 it 0.04, a quark mass ratio 

of ~ 21.2 lb 2.4, and a negligible departure from the Gell-Mann-Okubo mass 
[X^ '. formula. The mixing angles are 6*8 = -20.40° ±0.96° and Oq = -0.05° ±0.99°. 

We also give the values of all decay constants. It is shown that the common 
mixing pattern with one mixing angle Op is actually quite appropriate and 
algebraically related to the r]/r]' mixing pattern presently preferred by the 
ChPT community. For instance the traditional 9p is functionally related to 
the ChPT 9s and fulfills 9p ~ 6*3/2. The vanishing of supported by all 
data on radiative decays, gives a novel relation between mixing angles and the 
violation of nonet symmetry in the pseudoscalar sector. Finally, it is shown 
that the interplay of nonet symmetry breaking through U(3) — > SU(3)x U(l) 
satisfies all requirements of the physics of radiative decays without any need 
for additional glueballs. 
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I. INTRODUCTION 



We have recently proposed a model for radiative decays of all light mesons [|^] which gives 
a consistent and successful description of all reported experimental information. This covers 
the 14 decay modes of the kinclQy —>■ P7 and P —>■ '-ff. This vector meson dominance (VMD) 
based model relies on the hidden local symmetry (HLS) approach developed in Ref. 
which introduces the vector mesons as gauge bosons of a spontaneously broken hidden local 
symmetry and closely resembles Seiberg's EM duality in supersymmetric QCD as noted 
in Ref. and again in [^. Its anomalous sector (referred to hereafter as FKTUY), 
describes the radiative decays of light flavor mesons. In its original form this Lagrangian is 
U(3) symmetric, as it possesses both nonet symmetry and SU(3) flavor symmetry. 

In order to describe the full pattern of light mesons radiative decays, these unbroken 
schemes need to be supplemented with symmetry breaking mechanisms. Breaking the SU(3) 
flavor symmetry is an essential step [I]]. This is performed following the mechanism proposed 
by Bando, Kugo and Yamawaki (BKY) |^,^] and does not depend on any additional free 
parameter. An additional breaking procedure ||l|J^J^ is needed in order to describe the 
observed |T^ features of K* radiative decays; it allows us to recover for this sector a structure 



derived by G. Morpurgo in his approach to low energy QCD [pi]] . 

An explicit form |]T^ of nonet symmetry breaking (NSB) for vector (V) mesons seems 
to play a negligible role when focusing on radiative decays. Angular departures from 
ideal mixing are instead highly significant; they can be essentially explained by the uoj/ipi 
transitions^ inherent to SU(3) VMD models like HLS. One cannot, however, completely 
exclude that some kind of vector NSB is hidden inside these angular effects |P. 

Instead, NSB for pseudoscalar mesons (PS) is an essential ingredient [Q]. It has been 
performed in the manner of Ref. [T^, which turns out to allow couplings to singlet and octet 
PS components to be different. It is a purpose of the present paper to revisit the issue of 
how NSB can be consistently implemented within the HLS Lagrangian. 



The problem of rj/rj' mixing [23[l4| is tightly connected with the breaking of nonet 



symmetry. This can be performed at the level of the coupling constants [|r^, but this 



breaking can also be connected with a possible glue component inside light mesons |[T5|-[T8 
Indeed, the full set of radiative decays examplify the large importance of this effect. 
Within this context, this reference also showed that effects of such a glue component can 
only affect the rj' meson, but cannot be disentangled from genuine nonet symmetry breaking 
effects without some a priori knowledge of one of these twin phenomena. 

In connection with this particular problem, but even more closely related with the ef- 
fects of symmetry breakdown in Chiral Perturbation Theory (ChPT), Kaiser and Leutwyler 
||19| , |20[1 advocate anrj/rj' mixing scheme (see Ref. |^T|] for a comprehensive review), more com- 



plicated than the usual one, depending on two decay constants and two mixing angles. Some 



^ This counting does not include vr'' — > 77, which would serve to fix f-j^. We prefered using directly 
the PDG recommended value |10| as for fx- 

^ The subscript / indicates the ideal combinations. 
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phenomenological analyses p2|,^ have investigated this new scheme. In a more axiomatic 
approach to QCD, Shore |24] also finds appropriate a four parameter parametrization of the 
7]/ 7]' mixing. 

Nevertheless, the analysis of radiative decays (however, 14 independent decay modes) 
of Ref. does not find any need for a four-parameter structure of the rj/rj' mixing, as if 
phenomenology were exhibiting several relations among the 4 decay constants [|T^,ED|,E3 



which might be fulfilled at the (already high) level of accuracy permitted by the data. To 
be more precise, Ref. |T[] yields a quite satisfactory description of the data by introducing 
one PS mixing angle and one NSB parameter affecting the PS sector; additional departures 
from SU(3) flavor symmetry as per BKY arise only through a dependence in /x/Ztt, 
which can hardly be considered as a (free) parameter. 



More appealing, the anomalous Lagrangian of Wess, Zumino and Witten (WZW) p5| , p6 
with SU(3) symmetry broken as explained in Refs. |]T|, leads to deflnitions of the mixing 
angle and decay constants as per Current Algebra and as following from the HLS-FKTUY 
framework. The (single) mixing angle was found in Ref. to be ~ —10°, and, moreover, the 
value for the octet decay constant is /s = 0.82/^. The relevance of these parameter values 
is strongly supported by an impressive agreement with experimental data within a highly 
constrained model (5 parameters for 14 decay modes). Interestingly, Ref. ||2^, relying on 
lattice QCD calculations reaches also a mixing angle a value ~ —10° ± 2°, with a preference 
for -10.2°. 

On the other hand, there are repeated claims p8|-p0| that a (single) mixing angle, as 
coming from standard Current Algebra expressions, is quite appropriate and is found to 
be much less negative (~ —13° to ~ —15°) than expected from ChPT (~ —20°). A quite 
detailed discussion of this can be found in Ref. (see also |5^) where such an angle value 



is derived from a bound state approach. Another Lagrangian approach [Q, parent to HLS 
using a specific breaking scheme, recently claimed an angle value close to the previously 
mentioned ones (~ —15.4° ± 1.8°). 



All this seems in glaring disagreement with the expectations of ChPT ||r5| , p!4| , p!9| , |20| , |34| , |35 
Taking into account the special role of ChPT in low energy phenomenology, a possible 
contradiction between ChPT, lattice QCD calculations or the VMD conceptual framework]^ 
is a worrying question which must be addressed and understood. This is the purpose of the 
present paper, which will show that the contradiction is illusory and only due to different 
definitions of the same parameters in a naive understanding of the WZW approach (with 
encompasses the Current Algebra definitions) and in ChPT. We shall explicitly state the 
relationships between them. 

It will be shown that the VMD approach, relying on the HLS model broken as in Ref. [|l|, 
is actually in accord with all ChPT expectations associated with the ChPT Lagrangian 
||T9|^T1 This will be illustrated by deriving from a broken VMD Lagrangian model, 

all known leading order expressions for ChPT mixing angles and decay constants, and by 
deriving their expected numerical values. For the sake of conciseness, we shall frequently use 
NSB to refer to nonet symmetry breaking and to FSB for SU(3) fiavor symmetry breaking. 



'Or other approaches as listed in the throughout discussion in Ref. 
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The outline of the paper is as follows. In Section |T| we present a Lagrangian VMD 
model, based on the HLS approach, which includes both NSB and FSB. We show that there 
is a close connection between them. In Section |ITT|, we study the field transformation which 
permits us to write the kinetic energy of this broken VMD model in canonical form, in terms 
of renormalized fields. We show here that the field transformation of Refs. corresponds 
to a first order truncation in both NSB and FSB. 

Section |rV| gives the VMD description of the rj/r]' — > 77 decay which depends on one 
mixing angle and nonet symmetry breaking (as parametrized by x). In Section the 
corresponding description is derived, starting from the BKY broken WZW Lagrangian, and 
it is shown that WZW and VMD coincide. 

In Section 0, we derive the set of relations which allows one to define mixing angles and 
decay constants in accord with the standard (or extended) ChPT approach. Here we show, 
first that the definitions of mixing angles and decay constants from VMD /WZW and ChPT 
do not coincide once symmetry is broken and, second, that VMD provides expressions and 
values for all accessible ChPT parameters in accord with expectations. This is illustrated 
by several examples, including the functional relation between the VMD mixing angle 9p 
and the ChPT angle 6*3. 

In Section |V11| , we show that starting from the axial anomaly, it is possible to reconstruct 
the one angle mixing scheme as it arises in our broken VMD model and from the WZW 
Lagrangian; we comment on the previous use of ChPT predictions in phenomenological 
analyses of radiative decays data. 

In Section [VIII| , we show that nonet symmetry breaking and pseudoscalar mixing angle(s) 
are functionally related, which is a completely new result. This allows us to perform a fit 
of radiative decays with only 4 free parameters. The level of nonet symmetry breaking 
correlated with the fit value of the pseudoscalar mixing angle is shown to remove any need 
for glue in the 77' meson. A few other points of interest are also examined (quark mass ratio, 
isoscalar mass matrix, effects of NSB on PS mixing angle values). Finally, Section 11^ is 
devoted to conclusions. 



II. A BROKEN HLS MODEL FOR RADIATIVE DECAYS 

The model developed in Ref. in order to describe all light meson radiative decays 
relies on breaking nonet symmetry and flavor SU(3) in the HLS Lagrangian [Q], and especially 
in its anomalous (FKTUY) sector [Q. The breaking procedure of SU(3) flavor symmetry 
(referred to hereafter as FSB) in the non-anomalous HLS Lagrangian is the so-called new 
scheme, a variant the original BKY breaking mechanism [0] discussed in Ref. [Q. 

For the purpose of only studying light meson radiative decays [|l|, a detailed knowledge 
of the nonet symmetry breaking (NSB) mechanism is not needed; one only needs to know 
the field renormalization it would imply. The choice made in Ref. [|I| was to postulate a 
likely form; this was determined by the O'Donnell derivation of the SU(3) - not U(3) - 
VP'~^ couplings, which assumes only the SU(3) flavor group structure, gauge invariance and 



Lorentz invariance 12 



However, the way FSB and NSB in the PS sector merge together is a much stronger 
assumption which only relies on its impressive phenomenological success when describing 
the full set of radiative decays of light mesons. In this section, we aim at proposing a 
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Lagrangian model which provides the appropriate PS field renormalization; it allows to 
strongly motivate this assumption, by relating this Lagrangian to the ChPT framework. 



A. Basic Ingredients 

The basic ingredients of the effective Lagrangian approach to the interaction of vector 
and pseudoscalar mesons are the matrices V and P of the vector and pseudoscalar fields 
expressed in the fiavor {u, d, s) basis. The vector meson field matrix V is usually written^ 
in terms of ideally mixed states {u^ , 0^) 

p- (V + ^')/V2^*° , a = l,...,8. (1) 

K*- 



V = VT" = — 
V2 



V 



I 



Correspondingly, the pseudoscalar field matrix is usually defined as 



~ V2 



V 

8 



TT 



K- 



K+ 
3^8 + 73^0 



(2) 



using the conventional octet and singlet components (tts, r/o) for the isoscalar mesons. For 
definiteness, the SU(3) matrices will be denoted T" (a = 1, ■ ■ ■ 8) and fulfill the normalization 
condition Tr [T^'T^] = 6""^/ 2. We complete this matrix basis, by adding the unit matrix 
suitably normalized = 1/ 

The physical states {u, 0, rj, rj') are generated from the ideally mixed states by means 
of standard rotation angles 5v oi 5p for vector and pseudoscalar mesons. Correspondingly, 
the rotation angles for the singlet and octet states to the physically observed mesons are 
traditionally named Oy and 9p. These well known relations can be found in Refs. p|,p!2|,^,p!0 
The connection between ideal and physical uo and fields is treated heuristically in Ref. \ 
and rigorously in Ref. [0]. We recall for further use the traditional (one angle) expression 



(3) 



If fields undergo renormalization, the fields vrg and 770 in the expression above should be un- 
derstood renormalized . With a slightly liberalized, but obvious, notation, the expressions 
above for V and P can also be written 



' 7] ' 




cosOp 


— sin 6 p 






1 

.V . 




sinOp 


cosOp 




^0 



V = Vs + Vo, P = Ps + Po 



(4) 



which exhibit their octet and singlet component combinations, and show that nonet (U(3)) 
symmetry is implicitly assumed. 



The sign in front of 0^ means that we define 



\ss). 
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B. Physical Motivation for Nonet Symmetry Breaking (NSB) 



Referring to O'Donnell |T2 



V 



NSB implies modifying Eq. 

= V^ + yVo, P = Ps + xPo 



to 



(5) 



In this way, NSB changes the relative weight of the octet and singlet parts in a priori both 
meson sectors. Ref. has recently performed a throughout study of NSB in the vector 
sector and clearly concluded that data were consistent with no such NSB {i.e. ?/ = 1); it is 
therefore a motivated choice to neglect vector NSB and state y = 1 definitely. 

Another way to account for nonet symmetry breaking is to assume that the singlet sector 
contains a component other than the standard SU(3)/U(3) singlet; we name it glue only for 
convenience. A possible coupling of the rj/rj' doublet to glue can be accounted for by 
means of an additional angle 7, which is zero if one chooses to decouple this doublet from 
glue 







v' 




II 




.V 





COS 61 



— sin 6 p 







sin 6 



P cos 7 



cos6' 



P cos 7 sm 7 



— sin 6'p sin 7 — cos sin 7 cos 7 









Vo 




.99 _ 



(6) 



Indeed, following the analysis of Ref. |]T|, we do not have to introduce any coupling of 
the f] meson to glue, which would introduce an additional angle {j3 in Ref. The angle 7 
produces a coupling of (only) the rj' meson to glue. We have named rj" the possible triplet 
companion of the rj/r]' mesons, and do not attempt to identify it |. If 7 = 0, one clearly 
recovers the usual mixing pattern for the rj/rj' system by decoupling it from glue. 

When fitting the data on radiative decays of light mesons, the level of correlation between 
X and 7 is found such that assuming glue and exact nonet symmetry (x = 1), or assuming 
no glue (7 = 0) and some NSB (x ~ 0.9), provide the same description of the data 
Therefore, whether glue is required in order to describe the 77' properties is still a pending 
question which will be addressed in the present paper (see Section |V111D , when an educated 
guess about the value of x will be made. Let us note that the level of glue can be as large 
as 20% if nonet symmetry breaking is ignored []T],§ by setting x = 1 ; this conclusion has 



been reached also by others ||39|JT7|JT8[| . 

Clearly, in order that our conclusion on this point be of relevance, the exact meaning of 
X should be exhibited; a framework as reliable as ChPT is appropriate. This also motivates 
our goal of comparing our broken HLS-FKTUY framework to ChPT. 



C. Basics of the HLS Model 

We refer the reader to Ref. for a comprehensive review of the HLS model. A brief 
account can be found in Ref. IQ. We only recall the main features here. 



^We shall not also attempt to include this additional singlet in the Lagrangian model to be 
proposed for reasons which will become clear at the end of this paper. 
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The HLS Lagrangian can be written £hls = I^A + o-Cy-, where 



U = Tr 
P 



f2 

-^Tt[L-R]^ 
-yTr[L + i?]2 



(7) 



a is a parameter which is not fixed by the theory and f-,^ is the usual pion decay constant 
(92.41 MeV). The covariant derivative is 



(8) 



where is the photon field, the vector meson field matrix defined above and Q = 
Diag(2/3, —1/3, —1/3) is the quark charge matrix, e is the unit electric charge and g is the 
universal vector meson coupling 0. Finally, one generally chooses the "unitary" gauge, for 
which 



e« = a = e = exp(zP/A). 



(9) 



The standard VMD model is obtained by setting a = 2 in the HLS Lagrangian. However, 
several studies of the pion form factor |^,^ favor a ~ 2.4, quite inconsistent with 2. A 
simultaneous analysis of light meson radiative decays and vector meson leptonic decays [jl],|^ 
finds a ~ 2.4 — 2.5, quite consistent with pion form factor studies. 

The HLS Lagrangian is given in expanded form in Ref. (see Eq. (Al), where the 
pseudoscalar kinetic energy term has been omitted). For the purpose of the present paper, 
it should be noted that the pseudoscalar singlet field 7]q undergoes no interaction and only 
occurs in the (omitted) kinetic energy term. 



D. SU(3) Breaking Mechanism (FSB) of the HLS Model 

SU(3) symmetry breaking (FSB) of the HLS Lagrangian has been introduced by Bando, 
Kugo and Yamawaki |^ (already referred to as BKY) and originates from Refs. [^|^ . Brief 
accounts and some new developments can be found in Refs. connected more precisely 

with the anomalous sector . We refer the reader to Refs. [|T],^,|^,^,^ for detailed analyses of 
the properties of known variants of the BKY breaking scheme. Here we will only sketch the 
so-called new scheme detailed in Ref. [Q. Basically, the BKY breaking of SU(3) symmetry 
is performed by modifying Eqs. (^ in the following way 

J^A,v = -f Tr[(L tR){1 + (ae^.v-d + ^R^Ayd)m?- (10) 

which has a smooth unbroken limit. The constant matrices eA_y are given by Diag(0, 0, 
CAy)- Defining the breaking matrices X^ y = Diag(l, 1, 1 + CAy), these Lagrangian terms 
can be written 

CAy = -^Tr[(L t R)XAy{L T R)XAy]. (11) 
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The expanded expression of the BKY broken HLS Lagrangian can be found in Ref. 0] 
(see Eq. (A5) in the Appendix). One should note, among other properties of this breaking 
mechanism, that the pseudoscalar singlet field tjq does not undergo interactions with any of 
the other fields, as in the unbroken limit. It contributes only to the kinetic energy term in 

La = Tr [SPXa^PX^] + • • • (12) 

The basic consequence of this BKY breaking mechanism for FSB is thus to force a renor- 
malization of the (bare) pseudoscalar field matrix P, P ^ P' 

P' = X)l^PX]l\ (13) 



in order to restore the kinetic energy term to canonical form. The field transform in Eq. (|I3| 
has a smooth limit when Xa 1- Additionally, we have 




z = l + CA=\-^] = 1.495 ±0.030 , (14) 



The quantity z was named I a in Ref- ||1[- We shall also use the notation Z = l/2;~2/3 
in the following, for consistency with expressions written in Ref. |jl|. It should be noted 0, 
that the field renormalization (Eqs. (|TB|) and (P^) is required in order to recover the charge 
normalization condition F/^+(0) = 1. 

The correspondence with the usual ChPT Lagrangian is easy to establish. Indeed, con- 
cerning PS fields it turns out to consider the L5 term of the /Z'-^'* Lagrangian together with 
and the correspondence is z = 1 + SL^m\/ fl; this gives a quite good approximation 
of the expression for z = [fx/ f-wf and L5 ~ 2.14 10~^. 

Thus, the BKY breaking mechanism outlined above, results in a renormalization of the 
PS field matrix in clear correspondence with ChPT expectations. It does not result likewise 
in a renormalization of the vector field matrix [0,0. However, the correspondence between 



the results of Morpurgo |lll| and the so-called K* model ^ could well indicate that a 
renormalization of the vector fields is also needed 0, but only shows up in the K* sector. 
In the context of the present paper, we are actually independent of any kind of symmetry 
breaking in the vector sector. Let us only mention the main result of Ref. which tells that 
angular departures from ideal mixing is an appropriate parametrization of the system 
as long as one deals only with on-shell vector resonances. 



E. Nonet Symmetry Breaking (NSB) of the HLS Model 

We aim here at providing a reasonable mechanism for NSB in the PS sector within an 
effective Lagrangian. This was required in order to describe successfully the set of observed 
V^P7 and P77 radiative decays [Q. 

The main problem faced in the phenomenology of radiative decays is the generalization 
of Eq. (|13D to the case where NSB is also active. At leading order, this is solely determined 
by the influence of NSB on the kinetic energy part of an effective Lagrangian. 
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If FSB were absent, we already know that P ^ Pg + xPg is the required field renormal- 
ization, i.e. the renormalization results in a rescaling of the singlet part of the P matrix. 
This means that NSB should contribute specifically to the kinetic energy term which would 
become 

Ca = TT[dPdP] + cTildPodPo] ■ ■ ■ (15) 

in the absence of FSB. 

Let us now examine how we might incorporate the singlet contributions into the HLS 
Lagrangian. As is well known, the symmetry of the HLS Lagrangian is larger than 
SU(iV/) xSU(iV/), it is actually U(iV/) xU(A^/). However, this is unphysical. The extra vec- 
tor U(l) symmetry conserves baryon number and is thus desirable; moreover, as remarked 
above, this is supported by the data. However, the additional axial U(l) symmetry is a prob- 
lem as it would imply either parity doublets or a ninth light pseudoscalar (for reviews see 
Refs. [P^ -|45[| and recently I^Tj). Therefore, reducing the symmetry of the HLS Lagrangian 
is desirable. Introducing the chiral field U = C,iC,r = exp(i2P//^) 0, one obvious way is 
through determinant terms 

C = £hls + In det U ■ In det + In det d^U ■ In det d^U^ (16) 

where /i is a parameter with mass dimension and we have introduced the dimensionless 
parameter A to allow for nonet symmetry breaking. Considering the chiral transformation 
U — >• gl^UgR, we see Eq. (|16D is now only invariant under SU(A^/)xSU(A^/) or when ql = Qr 
(i.e., Uv), as desired. Rewriting the Lagrangian we have 

C = £jjLS + ^-^Tr In t/ ■ Tr In + A^Tr In d^U ■ Tr In d^Ul (17) 
Now recalling Eqs. (0) and (^, this can be rewritten 

C = Uls + C'^LS = -^HLS + ^fiVo + lxd^md^r]o (18) 

as Pq = ?7o1/a/6 and Tr[T^^^] = 0. Thus, through this breaking of the Ua{^) symmetry, 
the singlet acquires a mass which is nonvanishing in the chiral limit and an additional 
kinetic term. As can be clearly seen, this implementation of NSB only modifies the singlet 
contribution to the Lagrangian kinetic energy (and mass term) without changing the usual 
HLS interaction Lagrangian (see Eq. (Al) in Ref. [Q). 

It is quite interesting at this point to remark that the NSB parameter we introduce can 
be identified with the Ai coefficient of the C^^^ contribution to the ChPT Lagrangian |ll9| , |20 



as clear from Rel. (13) in who carries practically the same notations as ours. Therefore, 
the kinetic energy of the Lagrangian, which mostly determines the PS field renormalization, 
meets all expectations from ChPT. 

Having shown how U^(l) breaking might lead to an additional Lagrangian term, C as 
given in Eq. ([T^), we now wish to explore the consequences of this. We are interested in 
calculating the axial currents. This can be done through an infinitesimal (axial) variation 
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J^- - = f^Q^^ = V.^^[T^xxr'xAdP' + XfJ^'dr^^ (19) 

We see the octet components are unchanged, while the singlet component is affected by a 
factor of 1 + A. 



III. AN EFFECTIVE LAGRANGIAN MODEL WITH FSB AND NSB AT FIRST 

ORDER 

For the purpose of the present study, we are interested only in the PS kinetic energy part 



of the Lagrangian in Eq. (|T8D, which needs to be diagonalized in order to get the explicit 
transform P P', from bare to renormalized fields. Conversely, it is clear that, if NSB 
vanishes, the kinetic energy of the Lagrangian is rendered canonical by the transform in Eq. 
(O). Here, any reference to what can happen in the V sector is totally irrelevant. 



A. Diagonalization of the Effective Lagrangian Kinetic Energy 

The Lagrangian in Eq. (^) has a non-canonical kinetic energy, which is precisely of the 
form given in Eq. (|15|) with c = A. Putting it into a suitable diagonal form is thus required, 
in order to define the physical fields in terms of the unphysical (bare) field and get the axial 
currents in terms of the physical fields. 

It is suitable to perform diagonalization in two steps. The first step is simply to define 

1 /2 1 /2 

an intermediate renormalization step by P" = PXj( , which puts the nonet symmetric 
part of the kinetic energy term into canonical form. Practically, this means that pion fields 
are unchanged in this renormalization, while the kaon fields absorb a /x/Ztt factor, as if 
NSB were absent. Concerning isoscalar mesons, these (first step) renormalized fields can be 
expressed in terms of the original (bare) fields through 





= z 


B 


-A' 










-A 


C 




.^0 . 



The parameters A, B and C depend only on the FSB parameter z already defined and they 
are 

^=^(£^,0.16. B=(!i±i).0.90, C = <i±^ 0.80 , (21) 

3 z 2)Z 3z 

where the numerical values correspond to z ~ 3/2. A can be considered as the FSB charac- 
teristic size. C and B differ at first order in this breaking parameter since V^{B — C) = A. 
After this renormalization the kinetic energy T is still non-canonical. Using Eq. (pOD, T can 
be expressed in terms of the (intermediate, i. e. double prime) fields by 

2T = + [dri^]^ + X[Ad7i'^ + 59r/;f . (22) 

It is useful to define the FSB angle /3: 
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Diagonalizing Eq. (|22D gives the following renormalized fields 

vTg = COS /SvTg — sin [3rj'^ 

r]'o = [sin Pn'^ + cos pr]'^ ^1 + X{A^ + B^). (24) 

These field combinations, which directly follow from the eigensolutions of the quadratic 
form of Eq. (p2D, have a smooth limit when both FSB and NSB tend to zero (vTg vTg and 
Vo ~^ Vo)'i when FSB alone tends to zero {A — > 0) the limit is also smooth (vTg — > vrg and 
^0 ~^ ?7o\/l + A). However, the limit is not smooth when only NSB vanishes; indeed, Eq. ( ^if) 
shows that the fields remain rotated by an angle (5 which is non-zero if FSB is still active. In 
order to cure this disease, one can choose as final renormalized fields linear combinations of 
the solutions in Eqs. p^), which have the desired limit properties and conserve the canonical 
structure of T by the diagonalization. Using v = \Jl + \{A'^ + B'^) — 1, these combinations 
are 

VTg = (1 + f sin^ /9)vrg + v sin j3 cos Prj'^ 

% = V sin P cos /JvTg + (1 + t> cos^ P)r]Q (25) 



B. Field Transformation At First Order 

We can use directly the exact transformation given by Eqs. (p5|); however, our breaking 
procedure is actually leading order in all breaking parameters. This is illustrated by having 
shown the connection between our breaking procedure and the C^^^ + C^^^ ChPT Lagrangian. 
Therefore, it is meaningful to truncate the field tranform above at leading order. 

Truncating v at first order, we have v ~ X{A'^ + B'^)/2 ~ A5^/2, and then 



^0 



1 + -A' 
2 



A 



AB 



^-AB 
2 

^ 2 ' 























.Vo. 



(26) 



The last relation is obtained by removing breaking terms of order greater than 1. Using 
Eqs. (pop and (p6D, we can approximate the physical field (prime) combinations in terms of 
the bare fields by 







.V'o . 





B -A 



-A{1 + ^B') C{1 + ^B') 



Vo 



B 



-A 



Vo 



(27) 



This is the physical (first order) approximation which corresponds to the field renormal- 
ization used in Ref. and recalled in Eq. (pOj). The last matrix expression in Eq. ( pTf ) is 
obtained by remarking that XB^ differs from A by terms of order XA and then is legitimate 
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to neglect them at first order in further computations. The AX term in the lower leftmost 
matrix element is kept for consistency, but clearly plays a negligible role. 

Therefore, the field renormalization on which the study of Refs. relies is obtained 

from a Lagrangian model by truncating at first order in the breaking parameters. This 
provides an excellent fit to all light meson radiative decays, as can be seen from Ref. |I| 



and as will be shown below (see Subsection VIII B). We have checked that a fit to radiative 



decays performed as in Ref. [Q] using the exact field transformation in Eq. ( P^D instead of 
its first order approximation in Eq. gives indeed an improvement in fit quality, but 
negligible^. From these expressions, it is also clear that the NSB parameter x |T| is actually 

x = l- -B^ ~ , ^ ^ A ~ 0.20 - 0.25, (28) 
2 VI + A52 ^ ^ 



using the reference value for x (see Eq. (p^). We see the parameter A is small. In determining 
the accuracy and systematic errors, the neglected orders of magnitude should be estimated 
from the values of A [FSB] and (1 — x) [NSB]. It should be noted that, even if x carries 
prominently the information of NSB, it is somehow influenced by FSB as B = 1 + A/ \/2. 

In what follows we shall approximate the change of fields by its expression at first order 
in the breaking parameters, which can be written 

P = X^^/^(P^ + xP^)X^^/^ (29) 

The accuracy of this expression relative to the Lagrangian defined above can be estimated 
at ~ 5% by analyzing the magnitude of the neglected terms in Eq. 



IV. THE VMD DESCRIPTION OF rj/r]' 77 DECAYS 

Following FKTUY the anomalous U(3) symmetric Lagrangian describing PVV in- 
teractions is 

C = --^e'^'""'Tr[d^V,d,V^P]. (30) 

The PYj and P77 transitions amplitudes are obtained from this Lagrangian and the non- 
anomalous HLS Lagrangian, needed in order to describe the direct transition of vector 
mesons to photons. This non-anomalous HLS Lagrangian is given in its expanded form in 
Ref. [0 . It should only undergo the field renormalization of Eq. (^9]) , valid at first order in 
the (two) breaking parameters. 

The HLS model contains the Vector Meson Dominance (VMD) assumption (for a review 
see [^); it thus gives a way to relate the radiative decay modes VP'j to each other and to 
the P77 decays for light mesons, by giving a precise meaning to the equations sketched in 
Fig. I. 



^The improvement is larger if ones replaces the PDG value |10| for the rate (p ~^ v'l by the mean 
value of all presently available measurements ! 
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FIGURES 





V ,Y 




V 




FIG. 1. . 

Graphical representation of the relation among various kind of coupling constants. V 



and V stand for the lowest lying vector mesons {p^ 



); the internal vector meson lines 



are propagators at s 
squared. 



and are approximated by the corresponding tabulated [|T0| masses 



Propagating the field renormalization in Eq. 
Eq. (13) gives 



^ ) down to the FKTUY Lagrangian of 



-1/21 



(31) 



Then, the VVP Lagrangian is changed in a definite way by the renormalization procedure. 

The expanded form of the Lagrangian in Eq. (|3T|) is given in the Appendix of Ref . [|l| . 
The expressions for the two-photon decay widths of the rj and 1]' mesons can be derived 
from this and the non-anomalous Lagrangian. They are |]l| 



mi 



'-em 



cos t/p — v 2 X smUp 

3 3 

5-2Z . ^ + Z 

smt^P + v2 X cosUp 

3 3 



(32) 



'em 



77 



where Z = 1/ z = [f-w/ fx] . Actually, the last expression in Eq. ( p2[ ) is a normalization 
condition which allows us to fix the numerical coefficient in Eq. (|3T1) . It is clear that Eq. (|32D 
gives the two-photon decay widths in terms of /^r, fx-, x and only one mixing angle, 6p. This 
wiU be frequently referred to as the wave-function mixing angle (see Eq. (|^)). As normal 
these equations do not (and should not) depend on the vector meson parameters {g and 6v)- 
Thus, using standard Feynman rules, the HLS model provides definite expressions for 
the two-photon couplings of the pseudoscalar mesons, through its anomalous (FKTUY) 
sector. These expressions exhibit the traditional form [|1^,|1^,^] originally obtained through 
Current Algebra. These couplings are related to partial widths by 



r(X ^ 77) 



Ml 



X 



647r 



a 



1 



X 



vr 



(33) 



As a test, one can fit the parameters x and Op solely through radiative decays of the 
type VP'-) and use these values and their associated errors to predict the values for the two- 
photon decay widths of the 77 and rj' mesons. The fit values used for these computations ||l| 
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are x = 0.917±0.017 and 6p = —10.41° ± 1.21°. The results are given in Table 1 and clearly 
illustrate that the expressions in Eq. (^) are valid and that the VP'y processes accurately 
predict the two-photon decay widths. 

Stated otherwise, one does not need more than one angle {6p) in order to describe the rj 
and rj' radiative decays and this receives an especially strong support from all VP'j modes. 



Additionally, despite claims ||l3| , |35| , p6| , this angle is found to be ~ —10°, in apparent (as will 
be seen) inconsistency with the ChPT expectation of ^ —20° P^ , |55|J5^ . 



Mode 


VMD Fit 


PDG 


Comment 


Global Fit Quality 


{x,ep) 




Prediction 


Value 




X^/dof Probability 


Correlation 






0.514 ±0.026 


77 


11.07/10(35%) 


-0.34 


— s> 77 [keV] 


0.464 ±0.026 


0.46 ± 0.04 


PDG mean 


9.14/10(52%) 


-0.49 






0.324 ±0.046 


Primakoff 


14.82/10 (13%) 


-0.55 


T]' 77[keV] 


4.407 ±0.233 


4.27 ±0.19 


PDG mean 







Table 1 : Partial decay widths of the rj/rj' mesons, as reconstructed solely from fits to the 
radiative decays VP'y (leftmost data column) and their direct measurements [jlO| (second 
data column). The third data column displays fit quality parameters when using the 
corresponding rj measurement. The rightmost data column gives the correlation coefficient 
{x, 6p) in the corresponding case. 

A comment is of relevance concerning the data on 77 — » 77. One clearly sees that the 
VP7 modes considered altogether clearly prefer the PDG recommended (mean) value to ei- 
ther of the homogeneous reported measurements. Therefore, one may guess that the (single) 
Primakoff effect measurement and the (fourfold) 77 measurement, both suffer from system- 
atic errors in opposite directions. This guess is supported by the recent direct measurement 
of the ?7 — ^ 77 branching fraction ||^ 39.21% ± 0.3%, quite consistent with the PDG mean 
value. We shall revisit this issue in Section |VIII| . 

In order to substantiate the relative quality of the three data given in Table 1, we 
have redone the global fit, as described in Ref. |l[], changing only the 77 — > 77 data. The 
corresponding fit information is given in the rightmost pair of data columns in Table 1. Even 
if the fit probabilities are all quite acceptable, it is clear that the PDG recommanded value 
is indeed preferred by the full set of VP'-f decay modes. For this reason we use, from now 
on, the corresponding best fit results as reference values: 



X 



0.902 ±0.018, Oi 



-10.38° ±0.97° 



(34) 



V. THE WZW DESCRIPTION OF r?/^ ^ 77 DECAYS 



Starting from broken HLS and FKTUY, the VMD model of Ref. |]l| recovers the tra- 
ditional form for the two-photon decay amplitudes, {i.e. the one mixing angle expressions 
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of Current Algebra p^JT 
identification with our Eqs. 



Using these standard expressions, one indeed gets through 



k 

fs 



5-2Z 



k 

fo 



5 + Z 



6 



X 



(35) 



where Z = fxY, and /q § denote the (Current Algebra) singlet and octet decay constants; 
we have already defined 6p, the (single) mixing angle occurring in this approach. The /gg 
are named 77 decay constants in 

It is easy to check that Eqs. (0) and (|35| ) can be derived directly from the WZW 
Lagrangian [^,^. Indeed, this can be written 



WZW 



(36) 



(with Nc = 3) where Q = Diag(2/3,-l/3,-l/3) is the quark charge matrix, A is the electro- 
magnetic field and P is the bare pseudoscalar field matrix. Changing to the renormalized 
field P' through Eq. (|l3|) allows us to recover exactly the couplings in Eq. (|32D. 



This illustrates clearly that, what is named /§ in the Current Algebra expressions 
for rj/rj' decays to two photons, can be expressed solely in terms of and /i^-, in a way 
which fixes its value to /§ = 0.82/^. Correspondingly, we have /o = 1.17/^ which includes a 
correction of approximately 10% due to nonet symmetry breaking. The fact that the WZW 
Lagrangian leads to the same results as the FKTUY Lagrangian simply states their expected 
equivalence when deriving two-photon decay amplitudes. Stated otherwise, the structure of 
Eqs. ( P^ depends only on the BKY breaking Xa, with a small influence of PS NSB. 

However, the SU(3) sector of Chiral Perturbation Theory (ChPT) p^j3^j35[| is well known 
to predict fs/ fir — 1-25 and a mixing angle of ~ —20°. Then, the question is whether there is 
an inconsistency with respect to ChPT, or if there is a mismatch among definitions in ChPT 
(/o^s) and in the VMD/WZW approach (/qs); after symmetry breaking. More precisely, the 
question is whether /g and /q have actually the meaning of decay constants. 

Before closing this Section, it is interesting to compare our expressions for /q g with the 
corresponding ones in |pT| and [R^. Indeed, by identification of our Eqs. (B^) with Eqs. 



(26), (31) and (32) of Ref. |T[, we get 



T.o(0,0) 



fx 



2 
3 



(37) 



for the ratio of their reduced amplitudes. Their own numerical estimate for this ratio is 0.62, 
in quite good agreemeent with ours. In order to reach Eq. (pT]), we have stated x = 1 in our 
own expressions. So, one can consider that the present Eq. ( pS]) extends the results of Ref. 
II3TI (and ||3^) to the case when nonet symmetry is broken; the correction is however minor 
in this realm. 



VI. A CHPT DESCRIPTION OF r//r?' ^ 77 DECAYS 

We have seen above that, in the VMD procedure developed in Refs. the ex- 

pressions for /g, /o and the mixing angle are the same as those obtained from the matrix 
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elements for {'^'^\Cw zw\'u) ^^^d {'^'^ICwzwW) ■ We have denoted the parameters, obtained 
in this manner, /gg and 6p. In ChPT, however, the corresponding quantities (/o,8, ^8, ^o) 
are defined through other matrix elements, namely {^\d^ J^''^\ri) and {Q\d^ J^^\ri') , where the 
J*'° are the axial currents. It seems, however, traditionally admitted |ll3| , p5| , |5l| that both 



sets of definitions necessarily coincide. It is this last property which is addressed now. 



A. Usual ChPT Parameters from Broken VMD 

The axial current defined by Eq. (|19D, 

j; = -2U {Tr [r<^X^9^PX^] + d'^'^Xd.Po] } (38) 



can be rewritten in terms of the physical fields, through the transformation in Eq. (|29D (and 
also Eq. ( pTf ) for the isoscalar sector). We can write the matrix elements (0|J^|P'") for 
a = 1, ■ ■ ■ 7 and get the corresponding decay constants: 

{Q\j;"'\7^/K{q))=if^/Kq, (39) 

for pions and kaons, taking into account the expression for z. For the isoscalar sector, we 
get 



Jl = ^/.« + ^(1 - ^)/.«, 

Jl = ^(1 - ^)Ud,ni + ^/.x(l + A)9,<, (40) 

with obvious notations, and where we have used z = [fx/ J-kY' = 1/^- From above, we know 
that X = 1/VTT\B^, is infiuenced by FSB (B^ ~ 0.8). Moreover, the occurrence of simply 
(1 + A) - without any dependence upon z - is certainly due to the fact that NSB in the 
Lagrangian of Eq. ( |T8l) does not undergo SU(3) breaking effects. Therefore, it is consistent 
to consider that 1 + A ~ and make the (first order) approximation a;(l + A) = 1/x. 

One should note the occurrence in Eq. (|40|) of singlet field contributions to the octet axial 
current and, conversely, of octet field contribution into the singlet axial current. Additionally, 
these terms vanish in the limit of unbroken SU(3) fiavor symmetry {z = 1) as expected. 

These axial currents allow to define the following matrix elements 

(0J^|7r«(g)) = zfsq, , {0\J',\v\q)) = ^foQ, 

WMq)) = ^bsq, , (0| ^^^^(g)) = ^boq, (41) 

with 

/8 = ^^^^A = (1.33 ± 0.02)U , bs = y^{l- z)U = (-0.24 ± 0.01)/. 
fo = ^^^U = (1-29 ± 0.03)A , 6o = ^(1 - z)U = (-0.24 ± 0.01)/. (42) 
where the quoted errors are statistical only. So, at leading order in NSB, we have Bq = bg. 
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One readily observes a mismatch between the VMD/WZW definition for /§ and /o (see 
Eq. (|35|) ) and the ChPT definitions above; this mismatch is both algebraic and numerical. 
Otherwise, this /§ corresponds to the standard ChPT definition and it has its expected value 
l|,|3|,||,|5ll 



In order to switch to the matrix elements for {0\J'jl'^\ri/ri'), we use Eq. (^) together with 



the notations of Kaiser and Leutwyler |1T9| , [20[1 



\j'/\v/ri'{q)) 



0,8 



(43) 



and find 



V 



COS ^8 

F^ sin ^8 
-F° sin^o 
F° cos^o 



/8 cos^^p — ^8 sin 9p 
fs sin 6p + bs cos 9p 
bo cos 9p — /o sin 6'p 
bo sin 9p + fo cos 9p 



(1.269 ±0.008)/^ 
-(0.472 ±0.021)/^ 
(0.001 ±0.023)/^ 
(1.315 ±0.026)/^ 



(44) 



using the reference parameter values for x and 9p oi Eq. 



It should be stressed here 



that F'^/^ and 6'o/8 differ from /°/^ and 9p only by terms of order bo and b^. These cannot be 
neglected consistently if one keeps terms of order sin^p, which are numerically of the same 
order. Eqs. (PD lead to 



F^ = (1.36 ±0.01)/^ 
0s = -20.40° ± 0.96° 



F° = (1.32 ±0.03)/^ 
9o = -0.05° ± 0.99° 



(45) 



In the exact SU(3) limit {z = 1), the expressions in Eq. (|^) imply 9o = 9s = 9 p. Thus, 
the difference between them is, indeed, an effect of SU(3) flavor symmetry breaking, with 
only a marginal (numerical) influence of the nonet symmetry breaking parameter x. 

Now, the values given in Eq. (^) can indeed be compared with ChPT expectations, as 
these expressions correspond to the standard ChPT deflnition of mixing parameters. The 
value for F^ compares impressively to the parameter free prediction of Ref. pO| (1.34/7r with 
no quoted error). The value for F^ is harder to estimate theoretically because of its scale 
dependence; however, from the information given in Refs. |lT9| , pO| , F° ^ 1.3/7r seems in an 
acceptable range. The value ^8 = —20.40° ±0.96° is impressively consistent with all reported 
ChPT expectations (for example, -20° ± 4° from Ref. H, -20.5° from Ref. |19|). 

For ^0 the situation is unclear because the accuracy of the reported theoretical expec- 
tation [O — ~4° is lacking. Then, it is not possible to compare rigorously our result 



in Eq. ( ^51) with it; we show just below that the difference with ChPT (if any) is due to 
non-leading terms in breaking parameters. 



B. Further Comparison of VMD with ChPT 

One can ask about the correspondence of the expressions for the ChPT parameters 
coming from our VMD Lagrangian model of currents with their usual ChPT expressions in 
terms of /tt and fx- We show here that all expressions we can get in our HLS-FKTUY 
approach and all reported ChPT expectations to which they can be compared coincide 
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surely at leading order in breaking parameters. Let us illustrate with a few examples, 



mostly referred to in Ref. |119| , |21 
From Eqs. 



one can easily derive 



?Sl2 



[l + 2z) 



fl + 



2 <-2 



Let us define 



1 + 2e 



X 



1 + 5 



(46) 



(47) 



Neglecting terms of order e^, we have e = fx/ fn — ^ — 0.22 and 5 = x — 1 —0.10 as 
values for these perturbations to exact SU(3) flavor and nonet symmetries. Using Eq. 
Eq. ( ^61) gives 



[F 



812 



l + - + — + 0{e') fl 



which can be rewritten 



•ilF'f = ^fl - fl + 0{e\ 



(48) 



(49) 



This gives the same leading term as its ChPT expression (see Eq. (11) in Ref. |[T9[|). 
Concerning Eqs. (^) give at leading order 



[F 



012 



;i + y)(l-25) + 0(6^) 



which compares quite well with the Feldmann expression 

2/i + fl 



[F 



012 



from which one derives, neglecting terms of order 0{e5) 

K^ = -25= 4 - 1 = 0.20 ± 0.04 



(50) 



(51) 



(52) 



at the mass scale defined by radiative decays. One should note that this implies that Ai 
is nothing but what was named A as it should since nonet symmetry breaking is done in 
accordance with the ChPT Ci Lagrangian. Eq. (p2|) may serve to estimate, at the same 
scale, the other OZI violating parameters p!9| , pl| A2 and A3 from phenomenologyQiAs ~ 
-0.03 ± 0.02, A2 ~ 0.31 ± 0.02. These quantities are certainly scale dependent p|,|3, 
however this dependence is expected smooth . 
On the other hand, Eq. (HI) also gives 



' The symbol ~ is used because theoretical uncertainties of the relations among the Aj are not 
quoted. 
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F°F«sin {9s -9o) = fsbo + fobs- 
Using the expressions in Eq. ( ^21) and in Eq. (0), this is 



3F°F8sin(^^o-^8) = 2v^(/^-/^ 



(53) 



(54) 



which coincides at leading order with the corresponding quantity given in Ref. fl^ (see 
Eq. (13) there). One should note, however, that the leading correction is increased with 
respect to neglecting deviations from nonet symmetry, from 22% to 27%. On the other 
hand, one can check that 



2 



and 



4 



(55) 



in units of fj^, which also corresponds to the expectation that F° and differ only at first 
non-leading order Of course, leading and non-leading in our expressions refers to the 
small perturbation parameters of our model, e and 6 defined above. 

Therefore, at first order in breaking parameters, all expressions (and numerical values) 
deduced from our HLS broken Lagrangian meet all expectations of ChPT at the same order, 
namely, for F^, F^, 6s and 6q. Furthermore, using our fit [|| to radiative decays, we can 
provide Ai with a reliable value. 

One should however note that higher order corrections play some role; indeed, at leading 
order, F*^ and F^ differ by about 15% (see Eq. (^) just above), while the numerical value 
in Eq. (|45|) gives a difference of only ~ 3%. 

The single divergence - if any - concerns Oq. However, this is fully due to higher order 
corrections, the ChPT estimates of which are still missing. Indeed, this can be checked 
directly on Eq. ( p5[ ) which gives — 9s — —20°, while Eq. ( |5^ ) gives Oq — 9s — —15°, as 
obtained in Refs. |T9|j20[] , i.e. —14° to —16°. With this respect, it should be remarked that 
Eq. (Q) undergoes higher order corrections of more than 20%. 

The importance of having an improved theoretical estimate of 9^ should be stressed. 
Indeed, it should allow to reject several modellings (see Table 1 in Ref. |2l|), as the values 
proposed for 9o range between 0° and —30°. Our broken HLS-FKTUY framework points 
out that the condition — is well fulfilled by all data on radiative decays, as will be 
further illustrated in Subsection IVlllB. 



C. Relations between 6*8 , £^nd 9p 

Equations (|4^ ) and (^2]) together give 

/gtan^P + feg . fn , N . ^^(l-^) /c-RN 

taio.9s = - — — — = tan{9p + Lfp) , taiupp = V2 (56) 

fs — bstan9p (1 + 2^) 

which corresponds to ~ —10.02°, not infiuenced by NSB at this order. On the other hand, 
the same equations provide also 

tan 6^0 = y^""^f — ^ = tan {9p - ipp) , tam/jp = V2^j^-—^x (57) 
/o + oo tan 6'p (2 + -z) 
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Eqs. (|56|) and (^7^ together give 



es + eo = 2ep + ^(i - zf + ^(i -z){i- x) (58) 

up to terms of orders 0[{1 — z^] and 0[{1 — zY(l — x)], which represent a few percents only. 
This makes exphcit the connection between VMD/WZW and ChPT angles. Th. Feldmann 
has recently obtained a weaker form of the above relations in [21] (see its Eq. (84)), in the 
sense that NSB was neglected. 

Let us now consider the condition = as exact, which is clearly close to real life (see 
Eq. (^51) ). One can easily check that this condition can be cast into the form 

taxupp = KtanOp , (59) 

where K = (2 + z)/[(l + 2z)x] differs from unity by only 3%, which is likely well inside our 
(systematic) model errors. Therefore, Eq. (^) can be written 

tan6'8 = tan [9p + arctan (i^'tan6'p)] ~ tan26'p . (60) 

Thus we have 6p = 9^/2 already with an accuracy of order 1.5%. A slightly more accurate 
expression is obtained by expanding the relation above 

e^^{K + i)ep , (61) 

taking into account the (observed) smallness of Op. So, at the level of accuracy permitted 
by the data, there is a strict equivalence in using or Op. Furthermore, phenomenology 
indicates that higher order corrections should decrease the magnitude of 6*0. 



D. Estimate of the Quark Mass Ratio 

Since there is now some reason to trust the reliability of our numerical results in Eq. (^ 



we can use the following equation [19 



{[FX]' + [^.'M,f } = 4[/^M^]2^^ - [hM^n2S - 1) (62) 

in order to extract the ratio S of the strange quark mass to the mean value of the non-strange 
quark masses {S = rris/rh). This equation is actually second degree and thus potentially 
admits a spurious solution. The Mj terms are the corresponding meson masses. 

Using the mass values for neutral mesons and the information given in Ref. [|T^, we get 



— = 21.23 ± 2.42 or — = 2.51^-^ (63) 
m m 

The first solution compares well to the expectation of Current Algebra (25.9) and to the 



estimate (26.6) of Ref. [19| ; it is also in impressive agreement with the A. Pich estimate 
5B| , P7[ | 22.6 ± 3.3. The magnitude of the uncertainties (about 10%) is dominated by the 
errors on decay constants; errors due to choosing the neutral masses for the relevant mesons 
have not been accounted for. 
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VII. VMD, THE WZW LAGRANGIAN AND CHPT 



We have shown in Section |ITl| that the HLS model can be consistently extended in order 
to include nonet symmetry breaking along with SU(3) breaking effects. The connection 
with the ChPT Lagrangian £° + has been exhibited and proved successful. It has been 
shown reasonable to restrict to first order effects. The analysis of our HLS-FKTUY model 
with respect to ChPT, tells us that the BKY SU(3) breaking mechanism is justified. 
The most delicate assumption of the model for radiative decays of Ref. is the PS field 
renormalization in Eq. (pQ]); the diagonalization procedure above has shown that it is indeed 
perfectly justified at leading order in breaking parameters. 

Instead of radiative decays which exhibit a huge sensitivity |1|] to NSB (the parameter 
x), the basic ChPT parameters just considered exhibit only a marginal influence of NSB. 

Comparing the two sets of parameters, both derived from within a common VMD frame- 
work, the single clear conclusion is that there is a mismatch between the VMD/WZW cus- 
tomary definitions of decay constants and mixing angles and that currently stated within 
ChPT. Before commenting on phenomenological issues, we first show that this mismatch is 
a pure effect of SU(3) breaking which could have been foreseen since the work of Kaiser and 
Leutwyler [|l|,|g. 

A. From ChPT back to VMD/WZW 

The question of whether one can move back from the standard angles and decay 
constants of (extended) ChPT to the VMD/WZW framework is, of course, of special rel- 
evance. Indeed, we have already shown that, starting from our VMD/WZW model, the 
observed mixing angle of ~ —10° (for instance) was quite consistent will all expectations of 
ChPT, especially 9^ ~ —20°. The proof of the converse, i.e. deriving the WZW expressions 
from the axial anomaly, can be done on general grounds and is outlined in the Appendix. 
We detail here the algebra in order to illustrate the connection between ChPT concepts 
and the standard VMD parameters, and also test the consistency of having used first order 
approximation for the field transform. 

The basic idea is to remark that the divergence of axial currents is given by the axial 
anomaly at g — 0. In the case of two-photon decays, it takes the form 

(0|9^J^'l77) = iV(0|Tr[T'^F„^F"^]|77) , a = 0,3,8 (64) 

where the axial currents are given in Eq. (p^O]) , Fap is the photon field strength and Fap 
its dual, T°- are the SU(3) flavor matrices and is a normalization factor. Saturating the 
left-hand side of this expression with the nonet P of the (lightest) pseudoscalar mesons, we 
have 

(0|5,J^''^|77) = E(0|^M'^'''1^)7^G'p,, = iVTr[T'^g2] , a = 0, 3, 8 (65) 

where we have denoted by Gp^^ the decay amplitude of pseudoscalar mesons to two photons 
and Q is the quark charge matrix. 

As we limit the sum to the lowest pseudoscalar mesons, the single intermediate state for 
dJ^ is P = 7r° and then 
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(0|aj3|7r°) 



J-r 

77-u 



Tr [T"Q^^ 



(66) 



With (0|5J^|7r°) = /ttM^o, and because the last term in Eq. (|32D gives G.„a^^, Eq. (|66D 
provides the normahzation N = Ga/ir. 

For and 9J^, there are two possible intermediate states (the rj and t]' mesons) and 
Eq. (p5D gives: 



a 



(67) 



77 



Inverting these relations gives 
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-F^ + 2V2F^ 



(68) 



These expressions are the rj/r]' —>■ 77 amplitudes in terms of the Kaiser-Leutwyler pa- 
rameters and Oq/s- They can be reexpressed in terms of 9p and /o, /s, &05 by means 
of Eq. (HI) and become 
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G-q'j-y 



n^/S fofs — bobs 



{bo - 2V2/8) sin^p + (/o - 2V268) cos^p 
-(60 - 2V2/8) cos^p + (/o - 2v^68) sin^p 



(69) 



These expressions exhibit the standard Current Algebra structure. Having denoted the 
coefficients there by /085 we clearly have 



^ ^ fo - 2V2&8 
/s /o/s - bobs 



^ 1 feo-2v/2/8 

2V2 /o/s - 



(70) 



This proves that fs and fo would coincide with the rj/rj' decay constants /8 and fo , 
only if bo and ^8 were zero, i.e. if SU(3) were not broken ; in this case we would have 
fs = xfo = fn- Using the expressions in Eq. (|42|) and truncating at leading order in {z — 1) 
and (x — 1), it is easy to check that we get the fi as given in Eq. (|35D. 

The expressions in Eq. (|70|) are interesting in this regard: they clearly show that the 
mismatch originates from the fact that bo and ^8 are non-zero when SU(3) symmetry is 
broken {z ^ 1), which is basically the point of Refs. [|19|,^. Therefore, the usual expressions 
of Current Algebra do not directly give the isoscalar meson decay constants when SU(3) is 
broken. 

Nevertheless, this does not prevent Current Algebra from being the most economic for- 
mulation for the study of radiative decays, as it involves only two parameters {6p and x) 
instead of four highly correlated parameters (see Eqs. (^)). Moreover, we shall see in 
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the next Section how x and 9p are actually related by the observed smallness of 6q. This 
means that the Current Algebra formulation, suitably used, depends (practically) on a sin- 
gle parameter; this can be chosen equivalently as either of x (or Ai) or Op. In this case, 
the two-photon decay amplitudes for t) and rj' become a constrained system (2 equations, 1 
parameter) ! 

If one takes into account the smallness of expected from ChPT and from VMD esti- 
mates (see Eq. (0)), Eq. ( pSD provides quite an interesting result 

a 1 

sb's 

(71) 



G„-y-y + Gw^^tan^g — — j^-pz tt + 0{9q 

7rV3 -F8cos6'8 



together with Op = 08/2 + O(^o)- 



B. Phenomenological Consequences 

The consequences of the mismatch mentioned above are rather practical. Usually, the 
experimental determination of the pseudoscalar parameters (/o, fs, 6) is derived using the 
(quite standard) VMD/WZW Eqs. (0), which are nothing but the former equations of 
Current Algebra [|TB|,|T3|,|§ . 

However, generally, the model reconstruction quality ||38|,^ is defined by comparing fit 
results for /§, /o or Op with ChPT numerical expectations (F^, F*^, O^). This leads to highly 
confusing situations as shown by the discussions in Refs. [0,0 when justifying the value 
found for the mixing angle {Op ^ —13°) which mainly differs from our because of their 
neglecting NSB in the PS sector. The problem is strictly the same with decay constants 







Sometimes, numerical ChPT expectations for are attributed to what has been named 



/s in order to to constrain the WZW two-photon equations |ll3| , p!^ , |51| ; this mechanically 
(and artificially) pushes the mixing angle to ~ —20°, as can be seen from Fig. 1 in |T3 . 

What has been illustrated above is that such a phenomenological approach and such 
a theoretical treatment, are intrinsically inconsistent, as also noted by Ref. [^. Stated 
otherwise, it is meaningless to compare (or constrain) the Current Algebra equations (given 
also by Eqs. (0)) using the PS decay constants of ChPT or the value of the ChPT angle 
6^8, as traditionally done. 

To be more specific, we have proved that /§ = (0.82 ± 0.01)/,^, /o = (1-17 ± 0.02)/^ 
and a (single) Op ~ —10° which describes the rj/rj' mixing at the wave function level are 
consistently derived from VMD and/or the WZW Lagrangian after applying FSB (and 
NSB). Moreover, all this is perfectly consistent with all ChPT expectations at first order in 
breaking parameters: mainly ~ (1.25 — 1.35)/^ and ^8 — —20°. VMD has been able to 
provide new information (^o — and F^ ~ l-3/-n-, Ai ~ 0.20) of relevance for ChPT. A more 
refined comparison should wait until higher order ChPT estimates become available. 

As a side remark, one should also recall, from Ref. [Q, that SU(3) breaking does not 
affect the box anomalies for ^■n'^-n'rijrl . It is clear, however, from Eq. (39) in this reference, 
that nonet symmetry breaking can play some (numerically) minor role. Thus, all existing 
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analyses [^§,0 of the anomaly equations have to be redone from scratch, at least for 
consistency, knowing that the expected parameter values are not (directly) the ChPT ones. 
Related to this point, it is clear now that the Chanowitz equations correctly understood, 
do not point any longer to a failure of QCD, as incorrectly deduced in Ref. |^ because of 
the confusing angle problem mentioned above. 



VIII. FEEDBACK FROM THE CHPT PARAMETRIZATION 



The ChPT parameter values we have obtained (see Eq. (^)) allow for several remarks 
of importance which are to be discussed in this Section. 



A. A Hidden Relation between x and Bp /9s 



At the level of accuracy permitted by the whole set of radiative decays of light mesons, 
the results gathered in Eq. (^) indicate that = is well fulfilled experimentally. At its 
level of accuracy {Oq = —0.05° ± 1°), one can even ask oneself whether this relation is only 
approximate; this means that 9q does not undergo significant effects of SU(3) breaking, as 
opposed to ^8 and 6p. As remarked in Ref. [|l9l, this also means that, in the sense that |?7) 
is orthogonal to </°|0), the 77 meson is practically pure octet. But as shown above, the same 
\t]) happens also to be a mixture of |7r*) and |?7°) with an angle 6p ~ —10°. This illustrates 
the duality of definitions from another point of view. 

The numerical result ^0 = indicates that the state mixing angle 6p fulfills 

tan0p = v^|l^x (72) 
to good accuracy; this calls for several important remarks. 



By providing a definite value for x, Eq. ([7^) allows to address the issue of a possible 
glue content inside the rj' (see, for instance. Table II in Ref. or Fig. 2 in Ref. 
Indeed, this equation is purely a consequence of r] physics and Refs. |l|J^ have shown 
that no glue was required inside the rj meson. As a matter of consequence, Eq. ( [7^ ) is 
not influenced by a possible glue content in the t]' meson. 

Eq. ( |72D reveals an unexpected a/^eferaic relation^ between the mixing angles Op /6s and 
the nonet symmetry breaking parameter x. It will be checked explicitly in Subsection 



VIII B using the whole set of radiative decays. 



^ As a matter of consequences, there is a numerical correlation {x,6p) in fit procedures, which 
has been completely missed in our study in Ref. The correlation coefficient is given in Table 
1. This is smaller than what could have been expected, but the sign is consistent, if one remarks 
that 9p and x carry opposite signs. 
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Eq. (^) allows, for the first time, for a constrained fit to solely the 77 — 77 and 
rj — > 77 partial widths. Indeed, symmetry breaking effects in terms of singlet and octet 
components are completely determined by the BKY procedure, as shown in Section 
0. Only X is still free. Now, Eq. ([7^ ) tells that 6'p is not an independent additional 
parameter and is fixed by the value x. So, the single undetermined parameter is either 
Qp or X. 

We have performed the exercise and got Qp = —10.34° ± 0.22° (corresponding to 
X = 0.903 ± 0.017); using the PDG value of the r] rate - confirmed by the recent 
measurement of Ref. We get /dof = 0.8 10~^/1 (a 99% probability), while 



the Primakov measurement gives x^/dof = 6.99/1 (0.8% probability), and the 77 
measurement gives /dof = 3.89/1 (4.8% probability). So, the correlation is indeed 
observed; the reconstruction quality indicates however that the Primakoff measurement 
is affected by large systematics; this effect is present to a lesser extend in the mea- 
surement of the 77 Experiments. Therefore, a new accurate measurement of r/ — *• 77 
would be welcome. 

To our knowledge, it is the first time such a relation as Eq. ( [7^ ) is reported, li z = 1 
(then all 6''s vanish), x does not vanish but becomes unconstrained. Traditionally, the wave- 
function mixing angle is expressed in terms of PS meson masses ; our expression ( [72|) tells 
that the mixing angle can also be expressed as a function of fx/ f-w with some influence of 
NSB. 

From the point of view of ChPT, it should be stressed that Op is the variable which 
emphasizes the most the existence of nonet symmetry breaking; using the terminology of 
ChPT, it is directly proportional to x = 1/a/1 + Ai, which exhibits the scale dependence of 
Op. 



B. A New Global Fit to Radiative Decays 

From Eq. ([7^), we have been led to conclude that x and Op are algebraically related, at 
least to a very good approximation. In order to check this in its full realm, we have redone 
the fits given in Ref. requiring additionally this functional relation. This turns out to 
describe all radiative decays in terms of only 4 independent parameters {g, Op, Oy, (-t) or 
{g, X, Oy, £t), which is, by far, the most constrained fit of the 14 radiative decay modes 
ever attempted]^. We use the constant approximation for Sy following the conclusion of our 
study [|]. 

The fit quality obtained when setting up the constraint is /dof = 9.14/11, and does not 
exhibit any degradation compared to the fit quality reached when releasing this constraint 



now, we know that 9p and Og are functionally related, and that ~ —20° is equivalent to 
the favored |Q] 0p ~ —10°, we could, in principle, fix Og (and then Op) to its ChPT expectation. 
We have not perform this exercise, as the accuracy on the ChPT estimate of is still poor |34] 
and its sensitivity to NSB somewhat unclear. Nevertheless, it indicates that, from first principles, 
one can perform a fit to the 14 radiative decays, with remaining free fitting parameters referring 
only to vector mesons properties. 
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/dof = 9.13/10. The difference, in this last case, with Ref. [|l| is simply the use inside 
the fit of the PDG mean value for 77 — 77 instead of its mean value from all experiments 
(including the Primakoff effect measurement). In all cases, we have used the so-called |I| 
K* model, also commented on in Ref. 0. Therefore, Eq. (|72|) is indeed intrinsically present 
in the full data set examined sofar. 

Practically, the fit returns all parameters at the values obtained when leaving x and dp 
unrelated, even 9p which changes from Op = -10.38° ± 0.97° to Op = -10.32° ± 0.20°. The 
sharp reduction of the statistical error is an effect of removing the correlations by accounting 
explicitly for the functional relation in Eq. (|7^. The corresponding value for x is 0.901 (to 
be compared with the fit value x = 0.902 mentioned above). This numerical value has an 
important consequence, as will be seen in the next subsection. The fit branching fractions 
with X and 6p left free or related are given in Table 2 altogether with the data recommended 
by the PDG [|^, all used in the fit procedure. 

When setting the x — 9p relation, we do not observe any degradation in the quality of 
the description of the various branching fractions. Comparing the two sets of predictions in 
Table 2, one should note the sharp reduction of the statistical errors produced by having 
switched on the x — Op relation in all decay modes involving the r] meson. This also affects 
the modes involving though to a lesser extent. 

The increased accuracy we observe is not an artifact, but a trivial consequence of using 
independent (fit) variables instead of correlated ones. This explains why, even if it were 
theoretically better motivated, the use of correlated quantities like F^i^, and F^/r,' is not 
recommended in numerical analyses, as this results in unaccurate uncertainties. 

One should also note that recent measurements |Q of — >■ rj'"^ by the CMD-2 Collabo- 
ration at VEPP-2M, using new 77' decay modes, confirm the central value of Ref. [^] rather 
than that of Ref. reported by the PDG the agreement with the fit values we always 
get for this mode [l|] is thus improved. Indeed, the new measurement ([5.8 ±1.8] x 10~^) 
reported by Ref. (or by Ref. ||S^) would provide a minimum smaller than reported 
above by one unit; it has not been used in the fit in order to keep consistency with all 1998 



recommended values 10 



C. The x—9p Relation Kills Glue in the rj' Meson 

In the previous study of Refs. [|l],P it was shown that the correlation between glue and 
nonet symmetry breaking was huge. As stated several times above, accounting generally for 
glue coupling to the rj/r]' system implies that two angles have to be introduced in addition 
to 6p. One (/5) is such that /? = implies that the 77 meson does not couple to glue, the 
other (7) is such that 7 = implies a decoupling of the r]' from glue. 

However Table IV in Ref. [0 or Fig. 2 in Ref. clearly show that: i/ whatever the 
value of the nonet symmetry breaking parameter x, the angle (3 is not observed to deviate 
sensitively from zero; ii/ additionally, for x ~ 0.9 the angle 7 is also consistent with zero. 

We have seen above that it was indeed appropriate to parametrize PS NSB by x and given 
its most probable (fit) value 0.901. So, we can conclude that within the picture presented 
in this paper, there is no signal for a glue component inside the rj and the rj' mesons, but 
instead there is a significant signal for deviation from exact nonet symmetry: x = 0.901 
relative to 1 reveals a 5a significance level. 
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This is confirmed by performing tlie fit with 7 and x, now (numerically) decorrelated 
because of the functional relation in Eq. (^). In this case, the fit quality is strictly unchanged 
X^/dof= 9.14/10 and the minimum is reached for 7 = —0.02° ± 18°; this shows that no glue 
component inside the 77' is required by the data. 

The conclusion would be quite different p|,p|, p!7| , p^ , |39| if NSB could have been neglected. 
Indeed the level of correlation between standard nonet symmetry breaking (the parameter 
x) and glue is such that one can easily misidentify the former as being the later. So, even 
if our data set cannot exclude the existence of glue coupled to the rj/rj' system, it exclude 
presently its need. 



T>. 6p and the Isoscalar Mass Matrix 



In light of the above, the ChPT picture happens to be consistent with the standard one 
angle state mixing scheme of fields (or wave-functions): 
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The most accurate value for 6p comes out from fit to all radiative decays with the x — 9p 
correlation set up. Indeed, from the subsection just above, we know that it is legitimate to 
neglect coupling to glue. Actually, one cannot assert that glue (or any other singlet state) 
is not present inside the rj/rj' system, but what is shown in Refs. [|l],H for rj, and just above 
for rj', is that no glue contribution is required. A kind of minimum complexity argument, 
then leads to state /5 = 7 = and to decouple glue from the rj/rj' system. The angle value 
is 6p = —10.32° ± 0.20°, a hardly constrained value. 

With this, it is possible to revisit the determination of the isoscalar mass matrix 
11,11111 



Indeed we know that the mass matrix A4: 
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admits the following eigenvectors: 



Vr, = {cos 9p , — sin 9p) 
= (sin 6'p , cos 6'p) 



(75) 



with eigenvalues and M^,. This gives m 

^08 



0.320 ± 0.001, m 



2 

00 



0.898 ± 0.001 and 

m^g = -0.109 ± 0.004 in units of GeV^. This corresponds to moo"= 0.948 ± 0.001 GeV, 
and mgg = 0.566 ± 0.001 GeV, the former close to the rj' mass and the latter close to the 
rj mass as one might expect from the value of ^p. The off-diagonal term can be written 



m. 



08 



— 0.45M|^. The solution favored by VMD phenomenology is a very small deviation 



from the classical Gell-Mann-Okubo formula \u} 



)1|]. More precisely, from: 
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m;.= (1M^11 + A] (76) 

(with M^8 = russ) one extracts A ~ 0.01, where most of the error is due to choosing the 
mass values for K and vr. 



E. Broken VMD, ChPT and the Third Mixing Angle 

We have shown that the HLS model, after applying FSB and NSB, yields the structure 
of the ChPT Lagrangian C^^^ + C^^\ To be more specific, this concerns the PS kinetic energy 
part and the A2 and Lg terms [^,^ are not considered; they would influence the PS mass 



term which is outside the realm of this paper. Therefore, one can state that + 

is equivalent to the HLS Lagrangian broken as shown in Section |IT1| and the relationship 

exhibited in Section |VI C| between 6s and 6p is not accidental. 

Expressing the r]/r]' — >■ 77 coupling constants in terms of 6p, we have seen that this is the 
mixing angle which traditionally parametrizes the Current Algebra expressions. It is clearly 
a natural parametrization of VMD instead of the two mixing angles Oq and Og. For this 
purpose, it is quite interesting to compare our Eqs. (|32D with the corresponding Eqs.(39) in 
Ref. [^, expressed in terms of 9q and 6*8. 



Therefore, as it is possible to provide 6p a clear physical meaning through a mathe- 
matically well defined procedure, it is a quite legitimate choice. Its main virtue is to allow 
for a straigthforward handling of symmetry breaking effects in the VMD Lagrangians (HLS 
and FKTUY) as shown in Refs. 0,0]. The present paper has shown that the corresponding 
parametrization of VMD is easy to connect with the ChPT conceptual framework. 

F. Nonet Symmetry Breaking and the Third Mixing Angle 

Let us make a final (practical) remark on NSB. Ref. |Q, and Ref. |^ before, clearly 



proved that nonet symmetry breaking plays numerically a major role in accounting for 
radiative decays of light mesons within a relatively simple and constrained framework. Even 
if small in absolute magnitude {S ~ —0.10), the effect is statistically significant (~ 5a). At 
the present level of experimental accuracy, the data description quality is sensitive to this 
improvement. 

Quite interestingly, neglecting NSB in coupling expressions does not result in a dramatic 
change of the fitted mixing angle value, but mostly of the fit probabilities. Ref. [|l[ has thus 
obtained 9p = —14° ± 1° with, however, a degraded fit quality {x^/dof = 32/9). Table 1 in 
Ref. [m clearly show that all other physics parameters have unchanged fit values. 

This kind of angle value has been obtained in several other approaches ; their common 
feature is their neglecting NSB in the pseudoscalar sector. Ref. |2^ thus obtains 9p = 



-18.2° ± 1.4° in their fit to VP'y processes and 9p = -12.3° ± 2.0° in their fit to the 
P77 modes0. Refs. [pO| , pT| within the bound state approach and dealing with 2-photon 
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One should note that, even if both fit qualities are separately good, the two angle values are 
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processes only also get a mixing angle of 6p ^ —12°. The model of Ref. p3[ finds a fit 
solution 6'p ~ —15° (with no quoted fit quality) ; its freedom is however much larger than 
ours and the relation between their breaking scheme and ours - which now is motivated by 
ChPT concepts - is unclear. 

Instead, in our approach, the best solutions to separately VFj and P77 processes co- 
incide and both correspond to 6p = —10.3° ± 0.2°. Additionally, we should remind that, 
in our approach, all P77 couplings depend at leading order only on fx/ fir and 9p or the 
NSB parameter x. From this point of view, we consider that the result reported in (a 
prefered dp ~ —10.2°) gives our approach a strong support from lattice QCD estimates. 

It should thus be stressed that NSB in the pseudoscalar sector is exhibited, not so much 
by sharp changes in numerical values of physics parameters, but rather by their improved fit 
probabilities and, therefore, their actual accuracy. Indeed, as clear from this Section, and 
the previous ones, it was shown that the main ChPT parameters have values which are not 
sharply sensitive to having x 7^ 1. 

This is also the reason why the numerical results of Ref. |^ are so close of ours, even if 
these authors neglect NSB within the limited data set they consider. To be more specific, 
effects of the small value of 6 (or of Ai) are competing with SU(3) breaking, always by 
modifying non-leading corrections; for instance, the magnitude of the correction terms to 
9o — 9s (see Eq. (0), for instance). 



IX. CONCLUSION 

In a previous work, we were faced with a paradoxical problem. Using the HLS model 
and its anomalous FKTUY sector, together with a definite breaking scheme, it is possible 
to achieve quite a satisfactory description of all radiative decays, including 77/77' — >• 77. 
Within this framework, it was moreover possible to predict accurately these last rates, using 
only numerical information obtained by fitting the VP'y processes in isolation. This quite 
satisfactory pattern was obscured by some strange results: the (single) pseudoscalar mixing 
angle was found at 9p ~ —10° and the octet decay constant was /§ = 0.82/^, both in obvious 
disagreement with ChPT expectations. This meets independent analyses of other authors 
like Kekez, Klabucar, Scadron, Bramon ; it received recently a clear support of lattice QCD 
computations of the UKQCD group. 

The origin of this disagreement has been investigated. Starting from a broken VMD 
based Lagrangian, we have shown how to deduce the (WZW) rj/r]' two-photon amplitudes 
on the one hand, and the expectation values (0| J^'^l?]///') on the other hand (J's are the 
axial currents), which gives the customary (ChPT) definition of decay constants. 

The disagreement reported above has been traced back to inconsistent definitions for 
the same parameters provided by ChPT and the WZW Lagrangian (through the former 
definitions of Current Algebra). This inconsistency is a pure consequence of breaking SU(3) 
symmetry. For instance, it was shown that none of the two angles of extended ChPT can 



different enough that one can guess that the fit quality of a global description would be certainly 
degraded and providing 9p ~ —15°. 
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appear as such in the two-photon decay amphtudes. In some way, besides the angles 6s 
and 6q recently introduced by Kaiser and Leutwyler, the standard angle 6p, which still 
describes the 77/77' wave-function mixing, goes on playing an important role, the main one in 
radiative decays. It has been shown that 6p ~ 6^/2, within a ChPT motivated Lagrangian 
framework, perfectly accounting for all data on radiative decays and all accessible ChPT 
expectations. Moreover, analysis of numerical correlations has illustrated why the use of 
is more appropriate in order to get accurate measurements of physics parameters, including 
the standard ChPT ones. 

We thus have clearly illustrated that, the HLS-FKTUY model, supplemented with SU(3) 
flavor breaking a la BKY and nonet symmetry breaking, as was introduced in Ref. ^ was 
equivalent to the ChPT Lagrangian C^^^ + C^^\ Besides the successful description of all 
radiative decays, broken VMD thus meets all the requirements of extended ChPT (the two 
mixing angles and the two decay constants) with good accuracy. It was then shown that 
the relevant (and self-consistent) angle pattern is 6q ~ 0°, 6s — —20° and 6p ~ —10°. 
Correspondingly, we have simultaneously /§ = 0.82/,^ when using the WZW (or Current 
Algebra) definition and /§ = l.SS/jr, when using standard ChPT definition. Subsequently, 
we have shown that the cannot be interpreted as isoscalar meson decay constants because 
of flavor SU(3) breakdown. Therefore, VMD phenomenology is indeed able to provide ChPT 
with quite reliable input. With this respect, an important feedback from ChPT would be a 
more precise estimate of 60 and/or a (reliable) theoretical error. 

This study has led us to several additional conclusions: 

1/ Because ^0 = is observed with a quite impressive precision, it has been possible 
to relate nonet symmetry breaking (the parameter x) and the mixing angle 6p. To our 
knowledge such a relation has never been reported. 

ii/ The nonet symmetry breaking parameter A which weights the additional singlet 
contribution to the Lagrangian is small (A ~ 0.20 ± 0.04). It coincides with the usual OZI 
breaking parameter Ai of ChPT. 

ill/ As consequence of a; = 0.901, it has been shown, that no glue component is needed 
inside the rj' meson. 

iv/ By relating x and 6p, the condition 6q = allows to account for observed correlations 
in fitting radiative decays. Additionally, this leads us to propose a 4-parameter model to 
account for all data on radiative decays (including K*^ K^^)\ this is by far the most 
constrained model ever proposed and we proved that it is quite successful. 

v/ The quark mass ratio deduced from VMD information is rris/rh = 21.2 ± 2.4. 

vi/ Departure from the classical Gell-Mann-Okubo quadratic mass relation is observed 
at only the percent level. 
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Process 


X and 
Related 


X and 6p 
Unrelated 


PDG 


p 7r°7 (xlO^) 


5.16 ±0.03 


5.16 ±0.03 


6.8 ±1.7 


p 7r±7 (xlO^) 


5.12 ±0.03 


5.12 ±0.03 


4.5 ±0.5 


P ^ VI (xlO"^) 


3.16 ±0.05 


3.19 ±0.10 


o 4+0.8 


V ^P7 (xl02) 


33.3 ± 1.26 


34.5 ±2.1 


30.2 ± 1.3 


K*^ ^ K^-f (xlO^) 


9.80 ±0.94 


9.80 ±0.93 


9.9 ±0.9 


K*o ^ (xlO^) 


2.32 ±0.02 


2.32 ±0.02 


2.3 ±0.2 


u 7r°7 (xlO^) 


8.49 ±0.05 


8.49 ± 0.05 


8.5 ±0.5 


LO rj'y (xlO^) 


7.81 ±0.11 


7.88 ±0.23 


6.5 ±1.0 


1]' ^UJ-f (xl02) 


2.83 ±0.11 


2.94 ±0.19 


3.01 ±0.30 


7r°7 (xlO^) 


1.28 ±0.12 


1.27±0.12 


1.31 ±0.13 


0-^777 (xlO^) 


1.28 ±0.02 


1.27 ±0.04 


1.26 ±0.06 


(p^ri'-f (xlO^) 


0.59 ±0.02 


0.60 ±0.03 


1 O+0.7 


^ 77 (xlO^) 


38.87 ±0.75 


39.3 ± 1.8 


39.21 ±0.34 


r/' ^77 (xl02) 


2.09 ±0.08 


2.17±0.10 


2.11 ±0.13 



Table 2 : Radiative decay branching fractions. The first two data columns display the 
fit results using the K* model of Refs. [|l],§; in the first data column (present work) the 
x-6p of Eq. (|72|) is switched on while, in the second one, it is not. The data for rj j'-f 
is the recommended value |T^. The last data column displays the accepted values from 



the 1998 Review of Particle Properties |T0 . 
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APPENDIX A: THE ANOMALOUS DECAY TERM 



We shall now briefly discuss the equivalence of the anomalous decay amplitude T(P — *■ 
77) as calculated either from the divergence of the axial current, or the WZW Lagrangian 
p6|| . Let us recall the form of the axial current given in Eq.(|l9]). We then take the divergence 
of this to obtain]^ 

= 2fTT[T''XAT^XA]d^P^ + \f6''°d^P°. (Al) 
Now let us turn to the equations of motion for the pseudoscalar fields, namely, 

^ dC dC , . , 

leading to (allowing for a pseudoscalar mass term) 

2Tr[T"XAT^X^]a2p^ + A(5"°92p° = m^P" - ye^'^^^P^.P^^jTrfQ^T'^], (A3) 

where C is a well known dimensionless constant [^]. 

We are now in a position to show that the amplitude obtained form the axial cur- 
rent is equivalent to that obtained from the anomalous Lagrangian. First let us consider 
(AA|9^J^|0). As a total divergence this vanishes, in accordance with the fact that there 
are no truly massless particles in the spectrum (we are not considering the chiral limit). So 
using this along with Eqs. (^5) and ( |A3| ) we have 

/m^(AA|P"|0) = C{AA\F^M''Ti[Q^T'']\Q). (A4) 

As g — s> (X|P|0) = (X|P), hence, as rn?p is absorbed in the definition of the amplitude, 
we have 

T(P" ^ AA) = j{AA\F^^F^"'Tt[Q^T'']\0). (A5) 

We see that the result we have obtained, starting with the axial current, is the same as 
one would obtain from the anomalous Lagrangian term. This equivalence extends to VVP 
interactions. 



In this Section we use for conciseness the notation P = J2a=o —8 P°'T°'] thus for instance 
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